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ABSTRACT 

We present two-dimensional (2-D) radiative transfer modeling of the Eta Carinae 
binary system accounting for the presence of a wind-wind collision (WWC) cavity 
carved in the optically-thick wind of the primary star. By comparing synthetic line 
profiles with spectra obtained with the Hubble Space Teles cope /Space Telescope Imag- 
ing Spectrograph near apastron, we show that the WWC cavity has a strong influence 
on multi- wavelength diagnostics. This influence is regulated by the modification of the 
optical depth in the continuum and spectral lines. We find that Ha, H/3, and Fe II 
lines are the most affected by the WWC cavity, since they form over a large volume of 
the stellar wind of the primary. These spectral lines depend on latitude and azimuth 
since, according to the orientation of the cavity, different velocity regions of a spectral 
line are affected. For 2-D models with orientation corresponding to orbital inclination 
angle 110° < % < 140° and longitude of periastron 210° < oj < 330°, the blueshifted 
and zero- velocity regions of the line profiles are the most affected by the cavity. These 
orbital orientations are required to simultaneously fit the UV and optical spectrum of 
Eta Car around apastron, for a half-opening angle of the cavity in the range 50°-70°. 
We find that the excess P-Cygni absorption seen in Ha, H/3, and optical Fe 11 lines 
in I-D spherical models becomes much weaker or absent in the 2-D cavity models, in 
agreement with the observations. The observed UV spectrum of Eta Car is strongly 
dominated by absorption of Fe 11 lines that are superbly reproduced by our 2-D mod- 
els when the presence of the low-density WWC cavity is taken into account. Small 
discrepancies still remain, as the P-Cygni absorption of H7 and H<5 is overestimated 
by our 2-D models at apastron. We suggest that photoionization of the wind of the 
primary by the hot companion star is responsible for the weak absorption seen in these 
lines. Our CMFGEN models indicate that the primary star has a mass-loss rate of 
8.5 x 10~ 4 Mqyt" 1 and wind terminal velocity of 420 kins" 1 around the 2000-2001 
apastron. 

Key words: stars: atmospheres — stars: mass loss — stars: variables: other — 
supergiants — stars: individual (Eta Carinae) — stars: binaries 



1 INTRODUCTION 

The evolution and fate of massive stars is greatly affected 
by strong mass loss that oc curs via stellar winds through- 
out the star's lifetime (e . g.. lContilll976l ; lLanger et al.lll994l; 
iMevnet fe Maederl [2003). Rare giant eruptions of several 
solar masses may also occur during the unstable Lumi- 



* Based on observations made with the Hubble Space Telescope 
Imaging Spectrograph under programs 9420 and 9973. 



nous Blue Variable phase ( LBV; iHumphrevs fe Davidson! 
1 19941 : ISmith fe Owockil kooe). but the frequency and total 
mass lost via this phenomenon is uncertain. To add fur- 
ther uncertainty, the interaction with a binary companion 
has been suggested a s a possible trigger to giant eruption s 
ijKashi fe Sokerl 120101 : IPastorello etahl l201oi : ISmithl l201ll ). 
which highlights the concept that single and binary massive 
stars may have different mass loss histories and fates. Stud- 
ies of the most massive members of the LBV class, such as 
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Eta Carinae, are needed to gain insights into the mass loss 
processes in massive stars. 

Even though it is one of the most observed and 
scrutinized objects in the sky, the superluminous Eta 
Carinae still poses a great challenge to understanding the 
fate of massive stars. Located in the Trumpler 16 cluster 
in the Carina nebula and at a distan c e of 2 . 3 ±0 .1 kpc 
i|Walbornl 1 19731 : lAllen fc Hillierl 1 19931 : ISmithl 120061) , Eta 
Car is composed of a bright c entral object whose lu- 
mino si ty (L-,, ^ 5 x 10 6 Lq, iDavidson fe Humphreys! 
Il997l : ISmith et all 120031 1 ra nks among th e top in the 
Galaxy. First proposed by iDamineli et alj (|l997l ). it is 
now widely accepted that Eta Car is a massive, colliding- 
wind binary sy stem with a high eccentricity (e ~ 0.9, 
ICorcoranl 120050 and a n orbital period of 2022.7 ± 1.3 d 
|Damineli et alj|2008bl ). This scenario is supported by multi- 



wave l ength observations i n X- r avs (ICorcoran et al, 
200 ll. 120101: IPittard et"ail ll99Sl; I Pittard fc Corcorar 



Ishibashi et al 



2007 



1999; 



Henle y et al 



1997 



2002 



Corcoran 20051 : lHamaguchi et a* 
2008), ultraviolet (ISmith et al.l 



2004: llping et al l 120051: IMadura fc Grohl |2012|), optical 



Jvan Genderen et al. 2003 } 2006 : iFernandez-Laius et al 



20031 . [2009]. l20ld: ISteiner fc Daminelil |2004 iNielsen et al 



2007 • Damineli et alj l200d. l2008allbl: ITeodoro et all 



Mehner et al 



l201lt). near-infrared 
2004 iGull et alj 



Feast et al 



2012 



2001 
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Whitelock et al.l 12004 IGull et all 120081: . 
2010aHbl), and radio wavelengths (|Duncan fc White! 120031 : 
Abraham et al.ll2005l ). 

The orientation of the orbit has been the subject of 
debate over the last decade. Most authors agree on a lon- 
gitude of periastron of u ~ 24 0° — 270° , as proposed ini- 
tially by [Damineli et~afl ll 19971) . while others have found 



oj ~ 50° - 90° (e,. 



lAbraham et al] 120051 : iKashi fc Sokerl 
120081 ). Recently, IMadura et al l (|2012l ) determined the or- 
bital orientation for the first time in three-dimensional (3- 
D) space. Based on velocity- and spatially-resolved [Fe in] 
emission, these authors found that the Eta Car system has 
an orbital inclination angle i ~ 130° to 145°, longitude 
of periastron ui w 240° to 285°, and an orbital axis pro- 
jected on the sky at a position angle PA Z « 302° to 327° 
from North to East. These values are s upported by a multi- 
tude of multi-wavelength observatio ns (Okazaki et al. 2008; 
Parkin et al1l2009l. l201ll: iGroh et al.ll2010bl:lGull et alJboil ; 



Madura fc Grobll2012l ) 



The primary star (hereafter t/a) has been shown to be 
an LBV, with a high mass-loss rate of ~ 10 -3 Moyr -1 , 
wind terminal velocity of ~ 500 kms -1 , and dominates 
the luminosity of the system (|Hillier et alj I2OO1I . hereafter 
H01). Line profiles of Ha, H/3, and Fe n, obtained at dif- 
ferent positions in the Homunculus nebula that surrounds 
the binary system, support a latitude dependent wind only 
around apastron, with the wind beco ming more spherical 
during periastron |Smith et alj [20031 ). A scenario with a 
single, rapid-rotating 77A causing line profile variations as a 
function of latitude has been proposed to explain the deeper 
P-Cygni absorptio n profiles and fast er wind velocities ob- 
served at the pole |Smith et al.ll200"^ ). Rapid rotation may 
also explain the elongated K-band photos phere directly re- 
solved with interferometr ic measurements |van Boekel et all 
l2003l : IWeigelt~ et al.ll2007|). To explain th ese observations in a 
single-star scenario, IGroh et alj |2010al ) found that t]a must 
rotate between 77% to 92% of the critical velocity for break- 



up, with its rotational axis likely misaligned with the Ho- 
munculus polar axis. 

The nature and evolutionary state of the companion 
(hereafter ?7b), however, are far less certain since it has 
never been directly observed. This is because, in the visi- 
ble and longer wavelen gths, tja is a few orders of magni- 
tude brighter than ijb |Hillier et alj [20061 ). X-ray observa- 
tions sug gest a wind termina l velo ci ty of ~ 3000 k ms" 1 
for rm IPittard fc Corcoran! 120021 : lOkazaki et ail 120081 : 
IParkin et alj 20091 . l201ll ). which is in line with the detec- 
tion of high- velocity material (~ 2000 k ms -1 ) coming from 
the wind-wind col lision (WWC) zone jGroh et alj l2010bl : 
iHenlev et alj I2OO8T ) . Studies of the ionization effects of the 
companion on the nearby ejecta have yielded constraints on 
its temperature (T eB ~ 36, 000 - 41, 000 K iMehner et ail 
l2010l : IVerner et al1l2005l : ITeodoro et ail 12008!) and luminos- 
ity (10 5 £© < U < 10 6 £M. lMehner et alj|201ol ). 

The presence of a close, hot companion can potentially 
affect the integrated spectrum of Eta Car in two main ways. 
First, tjb, depending on its luminosity and temperature, 
may photoionize part of the wind of r]A, affecting H, He I, 
and Fe 11 lines, among others. This photoionization effect 
has been qualitatively discussed in various contexts (e. g 
Damineli etai1ll997l, |2008al : iHillier et alj|200il; [Smith etal 
2003; Davidson ct al. 2005; Nielsen et al. 2007; Gull et al 



20091 : IRichardson et ah l2010l : IMadura et alj|20ia ) 

Second, 3-D hydrodynamical simulations have shown 
that the wind of tjb greatly affects the geometry of the wind 
of t]a , creating a low-density WWC cavity in the wind of r/A 
and a thin, dense wind-wind interacting region between the 
two winds (IPittard fc Corcoran! 120021: Okazaki et al.1 [20081 : 
IParkin et alj|2009l . 1201 ll : IMadura et alj|2012l ). The effects of 
this WWC cavity on the spectrum of Eta Car have been 
qualitatively discussed in the literature, as a quantitative 
stu dy re quired yet unavailable 2-D radiative transfer mod- 
els. iHOll proposed that different lines will be affected differ- 
ently by rje depending o n the size of their formation region. 
IRichardson et ail (|2010| ) suggested that the carving of the 
wind of r/B might be responsible for explai ning variations 
in the Ha emission seen during periastron. ITeodoro et al.1 
(|2012h suggested that H<5 emission is not affected by the 
WWC during most of the orbit, while the Ha and Fe 11 
emissions are. 

Our goal in this paper is to perform a quantitative in- 
vestigation of the influence of r/B on the spectrum of tja, 
particularly due to the presence of a cavity in the wind of 
rjA- If these effects are significant, a detailed spectroscopic 
analysis might constrain the orbital parameters of the binary 
system and the individual stellar and wind parameters. Our 
ultimate goal is to gain insights on the nature of both stars. 
For instance, is t]a a rapid rotator, or can the presence of a 
close companion eq ually explain the available observations? 
IGroh et al.1 l|2010al ). based on the 2-D models discussed here, 
showed that the aforementioned if -band interferometric ob- 
servations of Eta Car can be explained by both scenarios. 

With these goals in mind, we apply, for the first time, 
2-D radiative transfer models to analyze multi-wavelength 
spectral regions of Eta Car, from the ultraviolet (UV) to the 
near-infrared. Here, we focus on the analysis of high spa- 
tial resolution spectroscopic observations obtained with the 
Hubble Space Telescope/Space Telescope Imaging Spectro- 
graph (HST/STIS) at UV and optical wavelengths around 



© 2012 RAS, MNRAS OOO,riJf20l 



On the influence of the companion star in Eta Carinae 3 



Table 1. HST/STIS Eta Car spectroscopic observations used in 
this paper 



Date 


0sp 


JD 


Wavelength Range 








(A) 


2000 Mar 23 


10.410 


2451626 


1170-1700 


2001 Apr 17 


10.603 


2452017 


1640-10100 



apastron. At this orbital phase, time-dependent and 3-D ef- 
fects due to orbital motion are minimized. 

We briefly describe in Sect.[2]the archival spectroscopic 
data of Eta Car used in this analysis. Section [3] presents a 
revised 1-D radiative transfer mo del of tja using the code 
CMFGEN ( jHillier fc Miller! 1 19981 ). Section H describes our 
2-D radiative transfer model for the wind of t/a including 
the WWC cavity and a dense WWC zone. In Section [5] we 
analyze how the presence of the cavity affects spectral lines 
of r/A during apastron. Section [6] presents a detailed com- 
parison between the 2-D model spectrum and the observed 
spectrum of the central source in Eta Car around apastron. 
We examine in Sect. [7| the influence of, and constraints to, 
the orbital parameters based on the analysis of the UV and 
optical spectra. The effects of the WWC cavity properties, 
such as the half-opening angle, are studied in Sect. [8] Effects 
of the WWC cavity on the continuum are discussed in Sect. 
[9] A discussion on the ionization effects expected from the 
hot companion is presented in Sect. \W\ and our concluding 
remarks are presented in Sect. 1111 

In forthcoming papers, we will investigate variations 
as a function of orbital phase, and the latitudinal and az- 
imuthal dependencies of the spectrum that can be probed 
using reflected spectra off the ejecta around Eta Car. 



2 OBSERVATIONS 

We employ UV and optical spectra of Eta Car compris- 
ing the inner ~ 0'.'2, minimizing the amount of contami- 
nation of the central source spectrum by the nearby ejecta. 
Seeing-limited ground-based observations are not adequate 
for our purposes as they are mainly contaminated by 1) 
bright, na rrow forbidden and permitted line emission from 
the eiecta JWeigelt fc Ebersbergeij[l98rj ; IWeigelt et alll 19951 ; 



the ejecta flWergelt fe IjJbe rsberger 198b; IWeigelt et 
iDavidson et alj Il995l ). and 2) faint, extended (^ 



0.3 



1.0") high-ionization forbidden line emission that ar i ses in 
the wind-wind inte racting region jGull et all [20091 . l201ll ; 
iMadura et 31112012! ). Both of these featu res, which are likely 
powered by ionizing photons from r/R JVerner et al.l 120051 ; 
iMehner et all 120101 ; IMadura et all 120121 ). are not modeled 
in this paper. 

The archival UV and optical HST/STIS observations 
analyzed in this paper were obtained by the HST Eta Car 
Treasury projeclQ A large wavelength coverage is needed to 
perform the detailed spectroscopic analysis, and reach the 
goals of this paper. Therefore, we use spectra gathered near 
apastron at <j) sp = 10.410 (2000 Mar 23) and <f) sp = 10.603 
(2001 Apr 17), where (j) sp corresponds to a given phase of 
the spectroscopic cycle of Eta Car. We adopt the ephemeris 



1 The reduced data are available online at http: / /ctacar. umn.edu 



from iDamineli et all d2008bl ) and the cycle labeling from 
iGroh fc Daminelil (|2004| ). 

Ultraviolet spectroscopic observations in the range 
1170-1700 A were obtained with the MAMA echelle grat- 
ing E140M, providing a spectral resolving power of R ~ 
48000 throughout this spectral region using an aperture of 
072 x 0" 2 Follow ing previous studies (|Nielsen et all 120051 ; 
IGroh et all l2010bl ). a 7-pixel (070875) extraction centered 
on the central source was adopted to minimize contamina- 
tion from the extended ejecta. Custom IDL routines devel- 
oped by one of us (JHG), assuming a boxcar function, were 
used to extract the spectrum. Optical and near-UV observa- 
tions were gathered with the 52" x 07l aperture, employing 
the G430M and G750M gratings to cover the spectral region 
1640-10100 A with R ~ 8000. An aperture of 6 half-pixels 
(07152) was used for the spectral ext raction with custom 
IDL routines. We refe r the reader t o Smith et all J2003T). 
IDavidson et~aH (120051) jHillier et all <|2006D , iNielsen et al l 
(|2007fi . and iGull et all (|2009T l for further details about the 
observations and data reduction. Table [T] summarizes the 
HST/STIS data used in this paper. 



3 THE 1-D SPHERICAL RADIATIVE 
TRANSFER MODEL OF rj A 

Since our 2-D models (Sectjffl are b ased on the spherically- 
symmetric models of 77A (|H0ll ; lH06t ) , we briefly summarize 
here the relevant model parameters obtained using the non- 
LTE, fully line blanket ed radiative transfer code CMFGEN 
(|Hillier fc Millerlll99SD . 

We use similar diagnostics and CMFGEN models as 

I — i 

those discussed by|H06J, with the key difference being a lower 
wind terminal velocity of Vao = 420 kms -1 and correspond- 
ingly lower mass-loss rate of M — 8.5 x 10 -4 MQyr -1 (Sect. 
13. 2|) . Because of the lower v^, several other parameters need 
to be adjusted to match the observed spectra. The parame- 
ters of our revised 1-D CMF GEN model are shown in Table 
[2] We note that, as in iHOll . models with R* in the range 
60 — 480 Rq produce similar fits to the observed spectrum. 
T he ato mic mode l and abundances are the same as described 
111 IHOll and H06, w ith the exceptio n that we adopt a solar 
Fe abundance from lAsplundl l|2005l ). 



3.1 Discrepancies when fitting the observed 
spectrum 

Although the spherical model of tja reproduces reasonably 
well the emission lines seen in the HST/S TIS s pectrum, 
some discrepancies remain (see Sect. 10 of IHOll) . Among 
these, the spherical model of Eta Car is not able to fit the UV 
and o ptical spectra simultaneously (Fig. 12 of iHillier et alj 
l200ll ). The revised CMFGEN models presented here have 
the same difficulties that we outline below (Fig. [1]). 

Significant improvements to the fitting of the UV spec- 
trum were obtained when including a "coron agraph" that 
block s the light of the inner 07033 region (|Hillier et al.l 
120061 ). In addition, the CMFGEN model also has difficul- 
ties with reproducing the central int ensity of both the high 
Balmer lines and Ha simultaneously (HOI). The optical He I 
lines are blueshifted compared to the 1-D CMFGEN model, 
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Figure 1. Spectrum observed with HST/STIS around apastron (black solid line) compared with the spherical CMFGEN model from 
iHOll with Voo = 500 kms -1 (blue dashed) and our revised spherical CMFGEN model with Voc = 420 kms -1 (red solid). 



which argues agains t them being formed in the wind of t)a 
|Nielsen et al.ll2007l ). 

The spherical model also overestimates the amount of 
absorption in the UV and optic al lines, in partic ular in H 
and Fe lines (Figs. 5, 6, and 7 of lHillier et~al1l200lh . This is 
surprising in a single star scenario, since CMFGEN mod- 
els are able to reproduce the H and Fe n P-Cygni pro- 
files of HDE 32682 5, which has an opti cal s pectrum simi- 
lar to that of ?7a (jHillier et all 1 19981 ). In lHOll . the spherical 
model of r]A was compared to observations obtained just 
after periastron (<f> Bp = 10.041). Since the amount of ob- 
served P-Cygni absorption is even lower ar o und apastron 
jDamineli et all 1 19981; iDavidson et al.l 120051 ; iNielsen et all 
120071 ; iDamineli et al.ll2008al l. the quality of the fits to the 
absorption profiles becomes much worse when comparing to 
the apastron data (Sect. [5]). The fits to the P-Cygni absorp- 
tion component and UV spectrum could be improved by 
significantly decreasing M in the models, whi ch conversely 
make s the fit to the emission lines much worse (|Hillier et al.l 
l200fj ). 



3.2 A lower wind terminal velocity for t)a 

When compared to HST/STIS observations obtained at 
apastron, we find that the original CMFGEN spherical 
model of r]A with Voo = 500 kms" 1 overestimates the 
amount of emission on the red side of most optical lines, such 
as H and Fe n lines (Fig. [TJ. In addition, that model over- 
estimates the maximum velocity of the P-Cygni absorption 
component of low-ionization resonance lines, such as C n 
A1334-1335, Si n A1526-1533, Mg n A2796-2802, and many 
resonance- like Fe n lines around A2400 and A2600 (Fig. [T]). 
We find that a spherical model with a wind terminal 



velocity of tja of = 420 kms" 1 provides improved fits 
to the UV and optical spectra. As shown in Figure [T] a 
model with this slightly lower value of Voo provides a much 
better fit to the red emission of optical lines of H and Fe II. 
Conversely, the fits to the blue side of the H emission become 
slightly worse. However, we prioritize the fitting of the red 
side of the lines profiles with the fore-knowledge that, for our 
preferred orbital orientation with i — 138° and u: = 260°, 
the red side of the line profiles are the least affected by the 
presence of the WWC cavity (Sect. [5]). The blue side of the 
optical line profiles are extremely affected not only by the 
WWC cavity as modeled in this paper, but perhaps also the 
driving of the wind of t?a might be modified by the ionizing 
photons of the hot, luminous t^b- 

We also analyzed the behavior of the maximum ve- 
locity of black absorption (hereafter Ubiack) of the low- 
ionization UV resonance lines of C 11, Si 11, Fe 11, and 
Mg 11, which has been demonst rated to be the best diag- 
nostic to Voo (|Prinia et al.l Il990t ) . Figure [T] shows that also 
the low-ionization UV resonance lines are more consistent 
with Voo ^ 420 kms" 1 than ~ 500 kms" 1 (Fig. [TJ). Our 
models suggest that, at least at apastron, material with ve- 
locities above ~ 420 kms" 1 appears only on the front side of 
the wind of t/a, and are probably ca used by turbulence b oth 
inherent to radiative-driven winds |Owocki et alj [l98&t ) or 
the WWC, with poss ible m ixing of the two winds. Similar 
to the 1-D models of IHOll . our models assume a microtur- 
bulent velocity of 30 kms" 1 that increases up to 50 kms" 1 
at larger distances in the wind of 17 a- 

The value of Voo obtained h ere is slightly lower than 
the Voo = 500 kms" 1 proposed bv lHOll based on HST/STIS 
data obtained just after periastron (ci S p = 10.041, 1998 
Mar). The change could be real or a result of using dif- 
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ferent diagnostics to determine Vao, and radiative transfer 
modeling near periastron is needed to properly address this 
issue. There are also issues related to the precise meaning 
of Woo, such as the influence of microturbulence and non- 
monotonicity arising from radiation line driving instabilities. 
At this point, given that Ubiack of the UV resonance lines 
does not change as a function of or b (|Hillier et al.l 120061 ; 
iGroh et al.1l2010bl ; iMehner et alj|201lf ). it seems more likely 
that u oo ~ 420 krns" 1 is more appropriate for material in 
the undisturbed wind of t/a at all phases. 

The value of Hoo found here based on detailed spec- 
troscopic analysis is lower than estimates using the same 
HST/ST1S data but ad-hoc assumptions without a quanti- 
tative physical mo del (vqq = 540 kms~ 1 JSmith et alj|2003t 
i-oo = 600 kms" 1 . IMehner et all 120111 ). The higher values 
of «oo from previous studies illustrate the challenge of ob- 
taining quantitative information about the velocity field of 
stellar winds without employing radiative transfer models. 

Interestingly, the modeling of high- ionization forbidden- 
line emission of 7/a, using an independent 3-D radiative 
transfer code and assuming Voo = 500 kms -1 , overestimated 
the amo unt of blueshifted e mission compared to the obser- 
vations l|Madura et al.|[2012f ). These authors suggested that 
a lower value of Uoo could possibly reconcile the 3-D mod- 
els and observations, and our results for the optical and UV 
lines of Eta Car provide further support for a lower Vao. 

Despite improving the fits to the line profiles in the 
UV and optical spectra, the discrepancies in the comparison 
with the observations (Sect. [3TTT> are still present. These are 
the main motivation for attempting a 2-D modeling of the 
spectrum of Eta Car, as we want to investigate whether the 
discrepancies found bv lHOll could be caused by the influence 
of tjb- 



4 THE 2-D RADIATIVE TRANSFER MODEL: 
PRIMARY WIND MODIFIED BY THE 
WIND- WIND COLLISION CAVITY 

We use an u p dated version of the 2-D code of 
iBusche fc Hillieij (|2005l . hereafter BH05) to create a low- 
density WWC cavity in the optically-thick, spherical wind 
of 7/a- Our 2-D model also includes the presence of a dense 
WWC zone, corresponding to the post-shocked primary 
wind. Figure [5] shows the as sumed geometry of the WWC 
cavity. We refer the r eader tolBH05l for further details on the 
orig inal code and to IGroh et al l (|200fl l200St 120091 , |2010al ) 
and iDriebe et all (|2009l ) for details on the code's applica- 
tions. Our implementation is flexible enough to be applied 
to other colliding-wind massive binaries, such as WR140. 

The 2-D code uses as input several quantities (e.g., 
energy-level populations, ionization structure, radiation 
field) from the spherically-symmetric model of 7/a computed 
with CMFGEN. The code allows for the specification of any 
arbitrary latitude-dependent variation of the wind density p 
and Voo. We initially consider a spherical wind for 7/a, with 
no intrinsic dependency of p and Voo on the colatitude 9. To 
create a low-density cavity in the wind of t/a, we modify the 
spherical p(9) so that the cavity has a conical surface with 
half-opening angle a (Fig. [2}. Material within the cavity is 
assigned a density b times lower than that of the spheri- 
cal wind model of 7/a. Thus, the parameter b reflects the 



different level populations and density of the wind of t/b, 
which fills in the cavity, with respect to the wind of t/a- At 
the moment, neither the flux from t/b nor the post-shocked 
secondary- wind are included. The WWC zone is approxi- 
mated by cone walls of angular thickness 8 a and a density 
contrast of f a times higher than the wind density given by 
the spherical mode l of t/a at a giv en radius. According to 
mass conservation |Gull et al.ll2009h . 

f a = [1 - cos(a)]/[sin(a)<y • (1) 

The simple conical shape is justified for orbital phases 
sufficiently away from periastron, when the ca vity has an ap- 
proximately 2-D axisymm etric conical form (jOkazaki et al.l 
|200S| ; iParkin et all [20091. Based on the expected location 
of the cone apex at a certain orbital phase and follow- 
ing the wind-wind momentum balance (|Canto et al.lll996l ; 
lOkazaki et al.l [2 0081 ), we place the apex of the cavity at a 
distance d ap cx from the primary star. The density struc- 
ture is axisymmetric about the x axis. Our implementation 
is flexible enough to use different parameterizations of the 
density st ructure as input. W e computed limited 2-D models 
using the ICanto et all l|l99rj ) analytic solution as input for 
the geometry of the cavity, which yielded similar results as 
those discussed in this paper. 

The 2-D models cover the same spatial scale as the 
1-D CMFGEN model, extending up to 5500 AU from rj A - 
Note that the long-term, large-scale effects caused by or- 
bital motion predict ed by 3-D hydrodynamical simulations 
(|Madura et all 120121 ') are not included in our models. To 
properly solve the radiative transfer equation and still be 
feasible to compute a grid of models, we use typically 70 
radial depth points and 51 angular points. Tests with up 
to 200 depth points and 90 angular points yielded similar 
results. 

The 2-D code of[§H05 allows the emissivities and opaci- 
ties to be either scaled or interpolated according to the mod- 
ified density. Both assumptions produce similar line profiles 
when including only the WWC cavity in the wind of Eta Car, 
and the precise choice does not affect the conclusions of this 
paper. Differences are only noticeable for the dense mate- 
rial in the WWC zone, in particular for f a > 4, when the 
interpolated source function predicts higher emission from 
optically-thick lines such as Ha. Hereafter, as a lower limit to 
the effects of the WWC zone, we present results for models 
with the scaled emissivities and opacities, except for the UV 
resonance and forbidden lines, where the interpolation of the 
source function is a more appropriate approximation. The 
observe r's frame spectrum is computed using CMF_FLUX 
|BH05l) . 

The geometry of the orbit is assumed as in 
iMadura et all (|2012l ). with t/b located in the x — y plane de- 
fined by the semi-major axis x and semi-minor axis y, and 
the orbital axis z being perpendicular to both. At apastron, 
t/b is located to the +x direction. The observer is located at 
an inclination angle i, which is the tilt angle between the or- 
bital axis and the line-of-sight (i = 0°: face-on view, i = 90° : 
edge-on view), and views the binary system according to a 
longitude of periastron u and at a certain orbital phase (/> or b ■ 
Periastron passage of the companion is defined as or b = 0.0 
and apastron as or b = 0.5. For simplicity, hereafter we as- 
sume that phase zero of the spectroscopic cycle coincides 
with periastron passage. Note that the precise timing of pe- 
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Figure 2. Left panel: Geometry of the Eta Car orbit (yellow ring in upper left inset) on the sky from lMadura et alj l|2012T ), The red arrow 
corresponds to the semi-major axis x, the green arrow to the semi-minor axis y, and the blue arrow to the orbital axis z. Middle and right 
panels: Assumed geometry of the wind of rj\ modified by the WWC cavity during apastron, showing a slice in the x — z plane (i.e., the 
binary system is viewed edge-on). The middle panel shows the large distance range covered by our 2-D radiative transfer computations 
(~ 5500 AU, ~ 2"), while the right panel shows a zoom into the inner regions of the WWC cavity. Both panels are color-coded to log 
density. 



riastron passage has not yet been constrained, but this does 
not affect our analysis of apastron observations. 

Our fiducial 2-D cavity model of Eta Car assumes the 
stellar and wind parameters of 7/a discussed in Sect.O while 
we adopt the orbital parameters and cavity orientation from 
iMadura et all (|2012h and references therein. We also assume 
the half-opening angle of the WWC cavity correspondent to 
a mo mentum ratio of r\ = 0.18 derived from the X-ray anal- 
ysis of lParkin et al.l (2011). Table [5] presents the parameters 
of our fiducial 2-D model that, unless noted otherwise, are 
assumed in the next Sections. 



5 EFFECTS OF THE WWC CAVITY ON 
SPECTRAL LINES OF VA 

Because the wind of Eta Car is extremely dense, the con- 
tinuum and line formation regions are much more extended 
than the size of the hydrostatic core (|H0l ; lH06l) . In addi- 
tion, the formation region of spectral lines is more extended 
than that of the continuum in the UV and optical. At these 
wavelengths, we anticipate that spectral lines will be more 
affected by the WWC cavity at apastron than the contin- 
uum. Here, we investigate thoroughly the reasons behind 
these effects since, in stellar winds, the formation of spec- 
tral lines is significantly more complex than the continuum 
formation. 

The formation region of spectral lines depends not only 
on the line opacity and emissivity, but on the energy levels 
involved in the transition, ionization potential (IP) of the 
species, non-LTE effects, and the line formation mechanism 
(i.e, recombination, resonance, continuum/line fluorescence, 
etc; see e.g., Mihalas 1978). In dense winds such as that from 
t)a, the ionization structure is stratified, and generally lines 
from species with high IP fo rm c loser to the star than lines 
from elements with low IP (|H0ll ). Therefore, spectral lines 
that are formed over a larger region of the wind of t/a are 



Table 2. Parameters of our fiducial 2-D cavity model of Eta Car 



?7a's Stellar and Wind parameters 


Value 


Ref. 


Luminosity (Lq) 


5 X 10 B 


DH97 


Temperature TV (K, at tr oss = 130) 


35200 


this work 


Temperature T c g (K, at tr oss =2/3) 


9400 


this work 


Radius R* (Rq, at rpt oss = 130) 


60 


H01 


Wind Terminal Velocity (km s _1 ) 


420 


this work 


Mass-Loss Rate (xl0~ 4 M yr _1 ) 


8.5 


this work 


Volume filling factor / 


0.1 


H01 


Distance to Eta Car d (kpc) 


2.3 


AH93, S06 


Cavity and orientation parameters 


Value 


Ref. 


Half-opening angle of the cavity a 


57° 


M12.P11 


Width of the cavity walls 8 a 


7° 


this work 


Overdensity of the cavity walls f a 


4 


this work 


Density ratio within the cavity b 


0.0016 


this work 


Distance d a pcx at apastron 


23 AU 


O08, M12 


Orbital Inclination i 


138° 


M12 


Longitude of Periastron u> 


260° 


M12 


P.A. on Sky of Orbital Axis PA 2 


312° 


M12 



1997), 



2008), 



Refer ences: DH 97=|Davidson fc Humphrey 

H01= lHillier et al.l feOOlD . O08=IOkazaJci*et al I 
M12= lMadura et alj <2012h . Pll= <Parkin et al.ll201ll) . 



expect ed to be more affected by the presence of the WWC 
cavity (|H0l[ ). Spectral lines formed in the innermost regions 
of 77a's wind, where the WWC cavity does not penetrate, 
will not be greatly affected. 

As an initial assessment of which spectral lines will be 
affected by the WWC cavity, we investigate the spatial re- 
gion where the flux from different spectral lines originate in 
the 1-D CMFGEN model of T]a- Recalling that, in a spheri- 
cally symmetric outflow, the flux in the observer's frame at 
a given frequency v is (e.g., Mihalas 1978) 
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Figure 3. Variation of the flux-like quantity pl(p) (color coded, normalized at each wavelength, and smoothed for clarity) as a function 
of impact parameter and wavelength for the 1-D CMFGEN model of va- The horizontal dashed lines in all panels correspond to the 
distance of the apex of the WWC cavity to rj\ for (from top to bottom) ipovb = 0.5, 0.9, 0.95, and 1.0. Spectral lines that have a 
significant fraction of pl(p) originating above the horizontal dashed line will likely be affected at that 4>orb- 



F »=tf J PUp)dp, (2) 

we investigate the behavior of the flux- like quantity pl v (p) 
as a function of impact parameter p, from the far-UV to 
the near-IR (Fig. [3|. This provides a first estimate of which 
lines will be affected by the cavity, since lines which have a 
significant fraction of their flux coming from p > d apex will 
be the ones more affected by the cavity. 

Figure [3] shows that the WWC cavity should affect the 



UV region even during apastron, especially in regions dom- 
inated by Fe n lines (e.g, AA1000 - 1250, AA1350 - 1390, 
AA1410 - 1460, AA2000 - 2600). This happens becaus e, for 
r?A, F e II lines are formed over very large spatial scales ()H0ll : 
HOG). In the optical and near-IR, the main lines that have 
formation regions larger than d apex at apastron are Ha, H/3, 
and Fe II lines. 

However, we stress that, because I(p) is a quantity pro- 
jected on the sky, the actual magnitude of the effects de- 
pends on the orientation of the WWC cavity. In addition, be- 
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v los 
[km/s] 



16.6 16.6 



228.8 383.0 



Figure 4. Line-of-sight velocity of the 1-D spherical CMFGEN 
model of t)a for an orbital orientation with i = 138° and lo = 
260°. Top: Three-dimensional rendering, with the thick black ar- 
row pointing toward the observer. The WWC cavity is shown in 
shades of gray. Note that this 3-D visualization does not corre- 
spond to the orientation on the sky, but it is rotated to better 
show the velocity regions of the wind of r}\ that are affected by 
the WWC cavity. Bottom: Cut in the x-z plane, with the yellow 
dashed arrow showing the direction to the observer, and the white 
solid line depicting the shape of the WWC cavity. 



cause the spectral lines are formed in an expanding outflow, 
different velocity regions of a spectral line form at different 
spatial regions (Fig. |3J). We find that spherical symmetry is 
broken by the presence of a WWC cavity and, depending on 
its orientation, different velocity regions of a spectral line 
are affected differently by the presence of low-density mate- 
rial within the cavity. As a consequence, the presence of the 
WWC cavity causes certain spectral lines to be latitude and 
azimuth dependent, in particular Ha, Hf3, and Fe 11 lines. 

Moreover, as t)b. moves in its eccentric orbit, d ape x and 
the orientation of the WWC cavity change. Therefore, spec- 
tral lines that are not expected to be affected by the WWC 
cavity at apastron, such as the higher Balmer and Paschen 
lines, could be affected at phases around periastron. This 
occurs because higher Balmer and Paschen lines form in 
the inner, undisturbed regions of the wind of n\ around 
apastron, where the WWC cavity does not penetrate. At 
later 4>orb, however, d avex is smaller and the WWC cavity 
likely disturbs the formation region of the higher Balmer 
and Paschen lines. In addition, because of the change in the 



WWC cavity orientation as a function of </f>ort>, lines that 
were affected at apastron at a given velocity interval will 
be affected at different velocity intervals around periastron. 
Ultimately, the presence of the WWC cavity will produce 
phase-locked variations in spectral lines. 

Here we focus on the effects of the WWC cavity at 
phases around apastron. For a WWC cavity orientation with 
i — 138° and ui = 260° , our 2-D models show that the wind 
of ?7b carves mainly the front-side of the wind of r]A, which 
contains material outflowing toward the observer (Fig. 3}. 
Only a relatively minor portion of ?ja's wind flowing tan- 
gentially to or away from the observer is affected by the 
WWC cavity. Therefore, the blue-shifted part of the line pro- 
files should be the most affected by the WWC cavity during 
apastron. 

Indeed, this scenario is confirmed by our 2-D radiative 
transfer calculations. For orbital orientations corresponding 
to the observer looking down the rarified cavity at apastron, 
the amount of Fe + in the line-of-sight is extremely reduced, 
causing less absorption of continuum radiation by Fe 11 tran- 
sitions. Therefore, the 2-D cavity model has much less iron 
blanketing than the 1-D CMFGEN model.This effect is most 
pronounced in spectral regions full of Fe 11 lines, such as 
around 1350-1450 A and 1800-2600 A. The extreme line 
overlap makes it complex to interpret individual features, 
which are usually blends of many spectral lines. 

Figure [S] presents monochromatic Ha images predicted 
by the 1-D and 2-D models at different velocities within 
the line. We investigate in detail the effects of the cavity 
for Ha, which is a recombination line formed over a large 
volume of the stellar wind. Analog effects are seen in the 
other unblended spectral lines, such as H/3 and Fe 11 lines. 

We find that the 1-D model has strong blueshifted ab- 
sorption at p smaller than the photospheric size at 6564 A 
and no emission from larger p (Fig. (5ji). In contrast, the 
2-D cavity model shows blueshifted emission at large p in 
directions toward where the cavity is projected on the sky 
(NW in this case, Fig. [5ji). A 2-D model including only the 
WWC cavity but no dense WWC zone (Fig. \Ejp) presents 
a similar behavior, showing that the key factor is the pres- 
ence of a WWC cavity in the wind of r?A- The enhanced 
blueshifted emission at large p is caused by the reduced op- 
tical depth because of the presence of the WWC cavity, al- 
lowing line photons that would otherwise be absorbed by the 
wind of t)a to escape. Since the flux at a give wavelength is 
F — J I(p)dpd5, the net result is a significant decrease in 
the amount of blueshifted Ha absorption. 

The WWC cavity also has a moderate effect in the emis- 
sion at zero velocities in Ha, allowing more line photons 
to escape toward the observer, but does not significantly 
change the image morphology (Figs. [5^,f). Because the back 
side of the primary wind is essentially not carved by the wind 
of r/B around apastron and remains undisturbed, no effects 
are seen in the redshifted emission of Ha when considering 
the WWC cavity alone (Fig.UJi). 

In addition to the WWC cavity, the line emission is also 
affected by the presence of a dense, shocked WWC zone. We 
find that, because of the orientation of the cavity, blueshifted 
Ha emission from the WWC zone is roughly aligned with the 
NE-SW direction, while redshifted emission is higher to the 
NW (Fig. [5f,i). This morphology resembles what has been 
seen in the extended forbidden line emission that originates 
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Figure 5. Monochromatic images of the Ho line predicted by the 1-D CMFGEN model (left column), 2-D cavity model with the 
parameters from Table [2] but without a dense WWC zone (i.e. f a = 1, middle), and the 2-D cavity model with the parameters from 
Table [2] The different rows correspond to representative images in the blueshifted (top row, v = —365 kms -1 ), zero velocity (middle 
row, v = kms -1 ), and redshifted (bottom row, v = 274 kms -1 ) parts of the Ha line profile. 



in the outer, time-averaged WWC zone (Gull et al. 2011), 
but scaled down by a factor of 30 spatially. 

Our models predict that Ha is optically thick, with a 
Sobolev optical depth larger than 1000 in the inner wind 
(for r < 200 AU). Because of that, the zero and redshifted 
images of Ha are also affected by absorption of line photons 
from the back side of the primary wind by the WWC zone. 
This effect is more pronounced for large f a . The absorption 
occurs mainly along the NE-SW direction (Fig.[S}i,i), and is 
caused by the redshifted line-of-sight velocity of material in 
the WWC zone toward these directions. 



6 COMPARISON WITH OBSERVATIONS 
AROUND APASTRON 

Building on the insights from previous Sections, we com- 
pare below the fiducial 2-D cavity model spectrum to the 
HST/STIS observations obtained at sp = 10.410 and 



10.603. As discussed below, our 2-D models show the pres- 
ence of a rarified WWC cavity in the dense wind of ija is 
sufficient for reproducing the UV and optical spectra of the 
central source, without the need of a rapid- rotating t)a- At 
this point, however, we cannot exclude that both scenarios 
occur simultaneously, with t]a being a rapid rotator. 



6.1 Ultraviolet spectrum 

The UV spectrum of Eta Car is strongly dominated by 
bound-bound t ransitions of Fe II believed to originate in 
the wind of tja (HOI; HOG). While both the r evised spherical 
model presented here and the one from lHOll have difficulties 
fitting the UV spectrum of Eta Car ( unles s light from the 
inner ~ 0'.'033 of the model is blocked; |H06|), our 2-D model 
reproduces reasonably well the UV spectrum (Fig. [5} . The 
quality of the fit improves compared to those from the spher- 
ical model because, for the orbital parameters assumed here 
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Figure 6. Comparison between the UV spectrum of Eta Car in the range 1260 - 3200 A observed with FST/STIS at <f> ap = 10.410 
(solid black line), the 2-D model including the cavity due to the wind-wind interaction (red dashed line), and the 1- D CMFGEN mode l 
of 7?A without the cavity (dotted blue line). Both models were reddened according to the interstellar extinction law of lFitzpatrickl (1999) 
with Ry = 4.0 and E(B — V) = 1.0, and scaled to match the observed continuum llHO ll) . 



(Table [5J , the observer looks down the low-density cavity, 
which dramatically reduces the column density, and thus 
opacity and optical depth of Fe II transitions in the line- 
of-sight to r/A- The end result is less absorption of Fe II in 
the UV spectrum while still maintaining the M of rjA at 
8.5 x 10~ 4 M0yr _1 . Therefore, according to our models, the 
UV spectrum is satisfactorily explained by the carving of 
the wind of rjA by r/B, without the need of a "c oronagraph" 
blocking the inner 0'.'033, as proposed by iHQq . Anomalous 
extinction could still occur in the line-of-sight to the central 
source (Sect. [TT]) . The fact that our 2-D models reproduce 
well the UV spectrum indicates that rjB has little flux contri- 
bution in the wavelength range analyzed here (1250-3300 A). 
This suggests that, perhaps, the luminosity of rjB is much 
lower than previously thought, maybe around a few times 
10 5 Lq a nd consistent with t he lower range of values sug- 
gested bv lMehner et alj ((2010) . Further updates of our code 



to include the flux contribution of t/b are needed to improve 
knowledge on rjB- 



Our model fails to reproduce high-ionization resonance 
lines, such as Si IV AA 1398, 1402 and C IV AA 1549, 1551 
(Fig. |6}. The P-Cygni absorption in t hese lines extends 
up to ~ —750 km s" 1 around apastron (|Hillier et all 1200(1 
iGroh et al.ll2010bl ). These velocities are consistent with the 
high-ionization resonance lines being formed in a region 
dominated by the turbulent velocity field of j^a's wind or 
in the WWC zone, perhaps moving sli ghtly faster t han th e 
terminal speed of rjA- As discussed in lGroh et all (|2010bl ), 
we speculate that these high-ionization resonance lines arise 
in modified regions of the wind of t/a, close to the WWC 
zone, that are photoionized by rjB (see Sect. I10|) . 
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Figure 7. Comparison of the spectrum of Eta Car observed with HST/ST1S in 2001 Apr 17 (</> sp = 10.603; black solid line) with the 
2-D CMFGEN models including the cavity due to the wind-wind interaction (red solid line) and the 1-D CMFGEN model of t]a without 
the cavity (dashed blue line). The fluxes are continuum- normalized. The main spectral lines are identified, with green and purple ticks 
marking the laboratory wavelength of H and Fe II transitions, respectively. 



6.2 Optical spectrum 



The assumed orbital and WWC cavity parameters could also 
be different (see Sects. [7] and |5J). 



The optical spectrum of the central ~ 0'.'2 of Eta Car 
is dominated by strong emission lines f rom H and Fe n 
|Davidson et al.lll995l ; Iffillier et al.ll200ll ) that were shown 
to ar ise in the wind of r/A ( HOll) . While the spherical model 
from lHOll strongly overestimates the amount of P-Cygni ab- 
sorption in these lines, our 2-D model yields much less P- 
Cygni absorption in Ha, H/3, FPy, and Fe n lines (Fig. 0). The 
explanation for the better fits is similar to that of the UV 
spectrum: since the observer looks down the cavity, there is 
less column density and optical depth in the line-of-sight. 
However, the 2-D synthetic spectrum still presents more ab- 
sorption than seen in the observed spectrum, in particular 
in H<5 and higher Balmer and Paschen lines. The absorp- 
tion part of these lines arises in regions of the primary wind 
unaffected by the WWC cavity at apastron, and thus 2-D 
and 1-D models have similar line profiles. A plausible ex- 
planation for the excess absorption is that the 2-D models 
currently neglect photoionization of the wind of t/a by t]b- 



Weak N II and He I lines are also present in the optical 
spectrum and most authors concur these are, at least dur- 
ing some orbita l phas e s, affected by photoionization of r/B 
(Damincli et al. 1997; Davidson et al. 2005; Nielsen et al. 
120071 ; iHumphrevs et all |2008| ; iMehner et aT1l201l[ ). In this 
case, the formation regions of these lines depend on the num- 
ber of ionizing photons from r\B, which depends mainly on 
its temperature, luminosity, and mass-loss rate. We found 
that the inclusion of the rarified WWC cavity and a dense 
WWC zone in the model is insufficient to provide better fits 
for these lines. As for the higher Balmer lines, the He I emis- 
sion and absorption that originates in the primary wind, if 
any, is unchanged by the presence of the WWC cavity alone 
at apastron (Fig. [7]). This is because the WWC cavity does 
not penetrate in the region where He I lines are formed. 
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7 CONSTRAINTS ON THE ORBITAL 
PARAMETERS 

As shown in Sect. O the UV and optical spectra of tja are 
severely affected by the presence of the rarified WWC cavity 
carved by j?b and the dense WWC zone. In this Section we 
analyze the constraints that can be imposed on the orbital 
inclination angle i and the longitude of periastron uj. For 
brevity, we quote results for selected spectral regions of ??a's 
spectrum: the near-UV (AA2400 - 2600), Ha, H/3, H5, and 
optical Fe n lines (AA4500 — 4600). Note that because we an- 
alyze spatially unresolved spectra from the central source, 
our models are insensitive to the value of PA Z . Interferomet- 
ric observations that spatially resolve s lines formed in th e 
stellar wind of t]a, such as Bracket 7 (|Weigelt et al.l 120071 ). 
could be used to further constrain PA Z . 



7.1 The orbital inclination angle i 

Our 2-D models indicate that the UV spectrum and the 
morphology of spectral lines significantly change as a func- 
tion of i. Figure [8] presents a set of 2-D models computed 
for different i, with the other parameters kept fixed (Ta- 
ble [2}. Orbital orientations with i and 180° — i produce 
identical model spectra, with the only difference being a 
counter-clockwise or clockwise orbital motion projected on 
the sky, respectively. Because the spectra analyzed here 
are spatially unresolved, this degeneracy cannot be broken 
with our models. Therefore, given the e vidence supporting 
a clockwise orb i tal motion on the s ky (|Smith et al.l |2004| ; 
iGull et alj|2009l ; iMadura et al.l [20121 ') . we investigate mod- 
els with 90° ^ i ^ 180°, noting that analogous constraints 
would be obtained for the parameter space 90° ^ i ^ 0° . 

As shown in Fig. [5] we find that most H and Fe 11 lines 
are latitude dependent. The strongest latitude-dependent ef- 
fects are seen in regions of the UV spectrum dominated by 
Fe 11 lines and in Ha, H/3, and optical Fe 11 lines. The differ- 
ences seen in the UV spectrum and strength of the P-Cygni 
absorption of optical lines as a function of i are both related 
to the varying amount of primary wind material in the line- 
of-sight, with the WWC cavity affecting different velocity 
regions of the line profiles (see Sect. [5} depending on i. 

The theoretical profiles presented in Fig. [^qualitatively 
resemble the observed line profiles scattered off different 
locatio ns in the Homunculus nebula shown in ISmith et alj 
(2003). Our 2-D cavity models predict that spectra scattered 
off the Homunculus poles should have stronger P-Cygni ab- 
sorption than those observed on the central star or equatorial 
positions. The presence of the WWC cavity thus mimic ob- 
served effects that could, in principle, also arise in a single, 
fast-rotating t]a- These have profound implications for in- 
terpreting the nature of the central source in Eta Car, and 
an in-depth investigation of this effect will be presented in 
an accompanying paper. 

Here we focus on constraining i, using as diagnostics 
the amount of Fe 11 absorption in the UV and blueshifted 
absorption of Ha, H/3, and Fe 11 lines. We find that orien- 
tations with 110° < i < 140° seem to be required to fit 
the UV spectrum and Ha line profiles satisfactorily (Fig. [5}. 
These values of allowed i correspond to orientations in which 
the observer looks down the low-density WWC cavity, with 
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Figure 8. Spectra from 2-D cavity models of Eta Car computed 
for different i, assuming a = 57° and lo = 260°. The panels 
show (a) the near-UV region around A2500, (b) H<5, (c) H/3, (d) 
Ha, and (e) optical Fe II lines in the region AA4500 — 4600. Note 
the latitudinal dependence on all spectral lines. In all panels, the 
black solid line corresponds to the observed HST/STIS spectrum 
at <4 SD = 10.603. 



the amount of Fe II and Ha P-Cygni absorptions reasonably 
matching the observations. 

The amount of primary wind material in the line-of- 
sight depends not only on i but also on uj and a. Variations 
in a could change the range of i that would fit the UV and 
optical apastron spectra according to 



110° < i < (90° + a), 



(3) 



for 25° < a < 90°. 

In addition to affecting the absorption profiles, chang- 
ing i might cause the emission line profiles to present sub- 
structures which vary in intensity and velocity. These are 
caused by the emission and absorption of dense, shocked 
material along the WWC zone (see Sect. [5] and Fig. [2J. 

7.2 Longitude of periastron uj 

We find that the presence of a WWC cavity induces not only 
latitudinal variations in some of the line profiles of Eta Car, 
but also azimuthal variations that can be used to constrain 
the longitude of periastron uj (Fig. [9} . Similar to the effects 
of i and a, the value of uj also regulates the amount of pri- 
mary wind material in the line-of-sight, and the amount of 
absorption in Ha, H/3, and Fe II lines, among others. 

Our 2-D models indicate that the best fit to the UV 
and optical spectra at apastron occurs for orientations with 
210° < uj < 330°, assuming a = 57° and i = 138°, corre- 
sponding to the observer viewing the binary system down 
the WWC cavity. These valu es are consiste n t with the 3-D 
orbital parameters found by IMadura et all (|2012l ). Chang- 
ing i and a could slightly affect the allowed ranges of uj, but 
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Figure 9. Similar to Fig- [HI but f° r different uj, assuming a = 57° 
and i = 138° . Note that the line profiles for d = 0° and 180° are 
very similar and are superimposed. 



will not allow models with . 
120051 ; iKashi fe Sokeill2008h . 

served spectrum. 



e.g. lj = 90° (|Abraham et al.l 
to be consistent with the ob- 



8 CONSTRAINTS ON PHYSICAL 

PROPERTIES OF THE WWC CAVITY 

In this Section we investigate how structural changes in the 
WWC cavity, such as the half-opening angle and density of 
the WWC zone, affect the emerging UV and optical spec- 
trum. In addition to obtaining insights on the physical prop- 
erties of the system, long-term changes in M and/or Voo of 
one or both stars might lead to changes in the structure of 
the WWC cavity, which in turn could affect the line profiles. 
For brevity, we present results for Ha, but similar effects are 
seen in other H and Fe n lines. 



8.1 The half-opening angle of the cavity a 

Figure 110a compares the HST/ST1S observations at (j) sp = 
10.410 with synthetic Ha profiles from our 2-D models for 
different values of a. The other model parameters are kept 
fixed (see Table 0. 

We find that changing a has a dramatic influence in the 
absorption and low-velocity emission parts of the P-Cygni 
profile of Ha. Similar to i and u, a also affects the amount 
of wind of r)\ that is carved by ?7b • The amount of carving 
determines whether the observer looks down the cavity or 
through the unmodified wind of rj a, dramatically impacting 
the Hydrogen n — 2 energy level population in the line-of- 
sight, which regulates the Ha absorption. We found that 2-D 
models with a < 50° show a strong, saturated P-Cygni ab- 
sorption component, since the observer's line-of-sight crosses 
mainly primary wind material (Fig. UOa V In contrast, for 
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Figure 10. (a): Ha line profiles computed with the 2-D cavity 
model assuming different a, for u) = 260° and i = 138°. The 
observed HST/STIS spectrum at rf> Bp = 10.603 is show in black. 
(b): Similar to (a), but for 2-D cavity models with different f a 
and <5 Q . 



models with a > 50°, the observer looks down the cavity 
and, consequently, Ha shows weak or no P-Cygni absorp- 
tion (Fig. 110b ). Also, as before, a change in i or oj would 
allow for slightly different a. Further constraints on a will 
be obtained in future studies of the variation of spectral lines 
as a function of <f> or b ■ 

The amount of Ha emission changes only slightly for 
models with a = 36° - 72° (Fig. ITPh). This effect occurs 
because the redshifted Ha emission comes essentially from 
the back side of the wind of r/Athat, for i — 138° and oj = 
260°, is only weakly carved by the WWC cavity (Fig. U). 
To see a significant effect in the emission strength, 77B would 
have to dominate the wind momentum balance in order to 
push the cavity to the back side of ?7a's wind (i.e., a > 90°). 

Our 2-D model favors 50° < a < 70° to fit the obser- 
vations since, as mentioned before, the observed Ha profile 
shows no P-Cygni absorption profile near apastron. A sim- 
ilar range of a is required to fit other optical lines and the 
UV spectrum simultaneously for tj — 260° and i — 138°. 
These values of a are consistent with those determined from 
analyses of the X-ray lightcurve ( Pittard fe Corcoran! |2002| ; 
lOkazaki et"al]|2008l ; [Parkin et al.ll2009T l. However, the value 
of M of r/A obtained here based on the UV and optical 
spectra (8.5 x 10 -4 Mq yr _1 ) is highe r than that derive d 
based on X-rays (4.8 x 10~ 4 M Q yr -1 . iParkin et al.l [201ll '). 
Consequently, our derived value of a requires rjB to have 
M ~ 2.0 x 10 -5 M0yr _1 in order to achieve similar momen- 
tum ratios as deduced from X-ray observations. 

In addition, the value of a could be biased due to ne- 
glecting the photoionization of the wind of r/A by t/b- In 
this scenario, portions of the wind of t]a which would have 
neutral H in the absence of 77B are actually photoionized, 
reducing the population of the n — 2 energy level of H and, 
consequently, the optical depth of Ha. This scenario might 
have important consequences, in particular if the observer's 
line of sight crosses a significant amount of the primary wind 
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close to the interaction region, which would be the case for 
Eta Car if indeed i ~ 138° and a ~ 57° . 



8.2 Density of the WWC zone f a 

Besides including microclumping via a volume filling fac- 
tor /, our 2-D models assume a WWC zone that is smooth 
on large scales (Fig. 0. Therefore, the calculations do not 
include large scale, optically-thick clumps, as would be ex- 
pected if the primary wind side of the shock cools radia- 
tivelv dPittard fe Corcoran! |2002l ; iParkin et all |2009| , 1201 ll ; 
iMoffat fe Corcoran! 20091 ). In addition, the wind of tja itself 
may contain large sc ale clumps, as suggested by polarimet- 
ric studies of LBVs l|Davies et all 120051 . [20071 '). Having that 
caveat in mind, we attempt here an exploratory investiga- 
tion of the spectral line morphology for different densities 
and widths of the primary side of the WWC zone. We as- 
sume that mass conservation ties f a , S a , and a according to 
Equation [T] 

Figure [TUb presents Ha line profiles illustrating the ef- 
fects of f a and S a on the spectrum of r/A- We note that 
similar effects are seen in other optically-thick lines affected 
by the WWC cavity, such as H,9 and strong Fe II lines. Our 
2-D models show that most of the effects of varying f a and 
5 a occur at low velocities, mostly on the redshifted side of 
the line profile. As f a increases, low-velocity redshifted ab- 
sorption from the back wall of the WWC zone affects more 
and more the Ha profile. Blueshifted emission and absorp- 
tion from the front wall of the WWC zone may also affect 
the Ha line profile. We find that models with f a ~ 1 to 4 
fit reasonably the Ha line profile at apastron, while model 
spectra with f a > 5 present sub-structures in the line profile 
that are not seen in the observed spectrum. 



40 



9 EFFECTS OF THE WWC CAVITY ON THE 
CONTINUUM OF rj A 

In this Section, we analyze how the presence of the WWC 
cavity affects the continuum of 77A during apastron. As we 
discuss below, these changes could occur due to t he 'bore- 
hole' phenomenon (Owocki & Smith, priv. comm.; [Madural 
2010) and modification of the optical depth of ?7a's wind, 
and em ission/absorption from the WWC zone (|Groh et al.l 
l2010al V 

As found bvlHOllandlHOr]. the continuum extension of 
r]A is regulated by the wavelength- and ra d ially- d ependent 
opacit y and emissivity . Madura fc Owocki (|2010T ) , iMadural 
|2010l ). iMadura et ail (^OllT ). and Madura et al. (2012, 
in prep.) showed that, depending on the opacity and 
emissivity of the wind of r]A, continuum flux changes 
at a given wavelength could occur when the apex of 
the WWC cavity reaches the surface of optical depth 
t = 2/3 at this wavelen gth. This 'bore-hole' phenomenon 
jMadura fc Owocki|[2oTol ) is likely related to the variations 
seen i n optical and near-IR photometry around periastro n 
(e. g., IWhitelock et al l 12004 iFernandez-Laius et al.l l2010h . 
We refer the reader to Madura et al. (2012, in prep.), where 
the continuum flux variation as a function of orbital phase 
is investigated in detail. 

Based on the revised 1-D CMFGEN model presented in 
Sect. 13.21 we investigate which wavelength regions could be 
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Figure 11. Distance from 77A where r = 2/3 in the continuum as 
a function of wavelength (black solid line). The horizontal dashed 
lines correspond to the distance of the apex of the WWC cavity to 
??A for </> rb = 0.5 (red), 0.9 (blue), 0.95 (brown), and 1.0 (green). 



affected by the cavity. Figure [TT] presents the distance from 
7]a where t — 2/3 as a function of wavelength for our revised 
spherical model. We find that the apex of the WWC cavity 
does not reach optically-thick regions of t/a's wind at most 
wavelengths for most of the orbital cycle. Figure [11] shows 
that at apastron the cavity reaches regions with r ^ 2/3 only 
for A < 912 A or A > 80 urn. Only at 0.95 < (t> rb < 1.05 does 
the cavity penetrate optically-thick regions in the optical 
and near infrared continuum. 

We now turn attention to the 2-D cavity models. Fig- 
ures I12b .b show the spectral energy distribution predicted 
by the 2-D cavity and the 1-D CMFGEN models from the 
far-UV to the near-IR. As anticipated from Fig. 111! our 2-D 
modeling reveals that the spectral energy distribution (SED) 
of t]a is strongly modified below 3000 A at apastron when 
the WWC cavity is open toward the observer (Figs. [T2b .bV 

Between 200-912 A, we find that a strong bore-hole 
effect may occur. The 2-D cavity model indicates that an 
observer viewing the binary system down the cavity toward 
7/a would detect a flux from 5 to 10 orders of magnitude 
higher than that arising from the 1-D model lacking a cav- 
ity. However, even with this extreme increase in flux, a hot, 
luminous companion star would easily dwarf the emission 
from the primary star between 200-912 A (Fig. [FZh). Be- 
cause of the steep radial dependence of the source function 
at these wavelengths, we anticipate that as the WWC cav- 
ity penetrates deeper and deeper at orbital phases closer to 
periastron, r/A will contribute an increasing amount of flux 
between 200-912 A. Whether this increase is enough to rival 
the flux emitted by r/B and produce significant photoioniza- 
tion of the ejecta depends on detailed radiative transfer cal- 
culations and on the properties of r]A, such as the assumed 
R*. 

From 3000 A up to 80 um, we find that the 2-D cavity 
model at apastron has an SED very similar to the one from 
the 1-D CMFGEN model, illustrating that no bore-hole ef- 
fect occurs around apastron in this wavelength range (Fig. 
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Figure 12. (a): Spectral energy distribution in the range 200— 
900 A of our fiducial 2-D cavity model of tta assuming the param- 
eters listed in Table [2] (red line) and the 1-D CMFGEN model of 
?7A (blue line). For comparison, the 1-D CMFGEN model of t)-q 
is also shown (green line), (b): Idem, for the wavelength range 
900-30000 A. All models were degraded to a spectral resolution 
of 1000 to help visualization in the UV and optical, (c): Flux ratio 
between the 2-D cavity model spectrum and the 1-D CMFGEN 
model spectrum. The red line corresponds to the flux ratio in- 
ferred for monochromatic continuum regions free of Fe II lines. 



112b) . This is consistent with the findings from, e. g, iMadural 
l|2010h and T. I. Madura et al. (2012, in prep.). 

Between 912-3000 A, however, there is a significant ex- 
tra amount of emission in the 2-D cavity model compared to 
the 1-D CMFGEN model fFig.H2bV even though the cavity 
does not penetrate into the optically-thick continuum re- 
gions of the wind of r)\ for these wavelengths (Fig. 1 1 1 p . We 
find that the increased continuum flux in the UV is caused by 
two effects: the modification of the continuum optical depth 
scale of ?7a's wind by the cavity, and the reduction in Fe II 
line blanketing. The latter effect clearly dominates the over- 
all flux increase, causing huge flux ratios of around 10-1000 
between 2-D and 1-D models within very narrow spectral re- 
gions (Fig. 112b). corresponding to (blends of) spectral lines. 
The effect of the WWC cavity on spectral lines is discussed 
in Sect[S] 

The increased continuum flux in the range 912- 
3000 A dilutes any putative emission lines from the compan- 
ion, making it challenging to directly detect spectral signa- 
tures from rjs- In particular, beyond ~ 1500 A, t)a strongly 
dominates the flux, which has important consequences for 
the il lumination of the circumbinary ejecta (|Madura fe Grohl 
l2012h . We also note that the observed spectral energy dis- 
tribution is different from the model predi ction s because of 
the anomalous reddening toward Eta Car (H01). 

We further investigate why around apastron there is an 
increase in flux at continuurrQ wavelengths between 912- 
3000 A and not at longer wavelengths, such as in the near- 
IR. Recalling that the intensity at a given frequency v and 
impact parameter p of a ray with source function S v (r) is 
(Mihalas 1978) 



(4) 



we study the behavior of r and S v (t) in the UV and near- 
IR. For simplicity, we consider a ray with impact parameter 
p = 0. 

Figure [T3k presents the variation of r in the UV (A = 
1384 A, red line) and near-IR (A = 21200 A, blue line) 
continua as a function of distance to tja for the 2-D cavity 
(dashed lines) and 1-D models (solid lines). We find that 
while the apex of the WWC cavity does not reach regions 
of the primary wind with r > 2/3 between 912 - 800000 A, 
the presence of the WWC cavity still modifies the continuum 
optical depth of ??a's wind at these wavelengths (Fig. 113a ). 
This occurs because material within the cavity has a much 
lower density, and thus much lower opacity Xi than that 
of the unmodified primary wind. Since in the usual p — z 
geometry (Mihalas 1978) 



-%dz , 



(5) 



the reduction in \ implies a reduction in the optical depth 
at a given distance z from the primary star. 

The variation of S u as a function of r in the UV and 
near-IR is shown in Fig. 113b . Our models indicate that S u (r) 
is steeper toward large r at UV wavelengths than in the near- 
IR. This property explains the larger relative variation in the 



I.e., regions not contaminated by Fe II lines. 
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Figure 13. (a): Radial variation of the continuum optical depth 
of a ray with p = (i.e., toward the core of tja) for the 1-D 
CMFGEN (solid lines) and 2-D cavity models (dashed lines). We 
show representative wavelengths in the UV (A1384, red lines) and 
near-IR (A21200, blue lines). Note that the dashed and solid lines 
are superimposed for r < 2 AU, as are the blue and red dashed 
lines for r > 8 AU. (b): Variation of the continuum source function 
of r) A in the UV (A1384) and near-IR (A21200) as a function of 
the respective continuum optical depth, (c): Idem, but for the 
intensity-like quantity Se~ T as a function of t. 



more noticeable in regions of high opacity /emissivity, such 
as in the UV (due to bound-free opacities from iron-group 
elements) and infrared (due to free- free processes). On one 
hand, the material in the WWC zone adds extra opacity 
in the line-of-sight to the continuum emitting region, which 
might decrease considerably the flux reaching the observer, 
depending on f a (Fig I14[) . On the other hand, the pres- 
ence of dense material in regions not covering the contin- 
uum source increases the emissivity of these regions. As a 
consequence, the amount of flux coming from large impact 
parameters will be higher. It is worth noting that emission 
and absorption from the WWC zone will be enhan ced at 
orbital phases around periastron (|Groh et al.ll2010al ). since 
the shock-compressed regions reach deeper into the wind of 

VA. 

An important prediction of our 2-D models, which in- 
clude the effects of 7]b via the presence of a WWC cavity and 
WWC zone, is that the morphology of the continuum and 
spectral line emitting regions should be dependent on the 
orbital phase. As can be seen in Fig. 1141 our models predict 
a roundish brightness distribution on the sky at apastron, 
with some asymmetry depending on the density and ioniza- 
tion structure of the walls of the WWC zone. Our models 
suggest that the image morphology should get increasingly 
more asym metric and elongat ed at orbital phases closer to 
periastron (|Groh et al.ll2010al ). This is in markedly contrast 
with the predicti on from a sin g le, rap id-rotator star scenario 
as proposed by ISmith et alj (|2003t ) . where rapid rotation 
durin g apastron should produce elongated, asymmetric im- 
ages (jGroh et al.ll2010al ). In addition, no time dependence in 
the image morphologies should be seen, unless the rotational 
velocity of t?a changes across the orbital cycle. Long-baseline 
optical interferometric observations are currently able to de- 
liver model-independent, reconstructed images that would 
allow one to test the aforementioned scenarios. 



intensity-like quantity S v {r)e~ T at p = in the UV com- 
pared to the near-IR (Fig. 113b ). Therefore, our 2-D models 
show that the change in the optical depth scale affects more 
the UV region than the optical/near-IR. Figure [T4l presents 
monochromatic images in the UV, optical, and near-IR con- 
tinua predicted by our 1-D and 2-D models, where the ef- 
fects of the modification of the optical depth scale can be 
visualized. 

Using similar reasoning, the high value of S v at large r 
in the UV also causes the 2-D cavity models to have higher 
values of I at large p compared to the 1-D model (Figs. 
114b . c). Ultimately, this results in a higher flux in the UV, 
of the order of 20-30% compared to the flux yielded by the 
spherical model (Figs. [12)) . In the optical/near-IR, however, 
little changes are seen (Figs. [T4"e .f.h.i). 

In addition, our 2-D model indicates that the pres- 
ence of the WWC cavity affects the intensity distribution 
of the images (and hence the flux) via emissio n from and 
absorpt ion by the shocked walls. As shown in iGroh et al.l 
(|2010al ) for the jf-band continuum, these effects can poten- 
tially change the morphology of the images and the amount 
of flux arising at large impact parameters. 

Emission and absorption from the WWC zone will be 



10 IONIZATION EFFECTS FROM THE 
COMPANION STAR 

At present, the 2-D code from lBH05l does not allow for the 
variation of the ionization stage and level populations as a 
function of latitude. Therefore, material inside the WWC 
cavity and along its shocked walls have the same ionization 
structure as the wind of 77a- This is admittedly a zero-order 
approximation in the case of the Eta Car binary system, 
mainly because of the two factors elaborated below. 

First, r/B supposedly has T e s = 36,000 - 41,000 K 
l|Verner et al.ll2005l ; iTeodoro et ai1l2008l ; iMehner et al1l2O10D 
and, depending on its luminosity, might be able to ionize hy- 
drogen and other species in a non-negligible fraction of the 
wind of tja ■ This could affect the P-Cygni absorption of hy- 
drogen lines, and give rise to spectral lines of high-ionization 
stages (such as Si ivAA 1398, 1402) that would otherwise not 
be present in the sp ectrum of tja in a single star scenario. 
iMadura et alj (|2012l ) estimated the volume of the wind of t)a 
that is photoionized by t/b and found that a significant frac- 
tion of H could be ionized near a pastron (see t heir Fig. 6). 
This result has been confirmed by iKruipl (|201ll ) using inde- 
pendent radiative transfer modeling. However, as a caveat, 
both calculations include only H and He, while iHOll showed 
that line blanketing by the Fe group elements is a key factor 
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Figure 14. Monochromatic continuum images predicted by the 1-D CMFGEN model (left column), 2-D cavity model with the parameters 
from Table [2] (middle) , and 2-D cavity model with a denser WWC zone with f a = 10 (left). The top row shows the UV region at A1384, 
the optical (A5500) is presented in the middle row, while the near-IR region (A21200) is displayed in the bottom row. 



to the H ionization structure in dense stellar winds, as is the 
case of rjA- 

Second, material is shocked along the wall, significantly 
increasing its temperature, and is compressed and radia- 
tively cools on the side of the primary. In addition, ionizing 
photons from r]B illuminate the post-shocked primary wind, 
and will affect this material. Hot, post-shocked secondary 
wind material will also be i n line - of-sigh t and likely emit- 
ting in X-rays iParkin et ail (|2009l . |201lh . Because of these 
complications we are unable to self-consistently model the 
line emission and absorption from the walls of the shock 
cone. However, they do not affect the conclusions reached in 
this paper regarding the WWC cavity effects. We note that 
the effects described above should affect much less the con- 
tinuum emission and absorp t ion, in particular the near-IR 
region studied in lGroh et al.l (|2010al ). 

As a consequence of these limitations, our model is not 
able to reproduce high-ionization lines thought to arise from 
material photoionized by r/B- We defer a detailed model of 
ionization effects to future work, when it will become feasi- 
ble to couple 3-D radiative transfer tools (e.g., SIMPLEX; 



Paardcko oper et al. 1 120101 ; lKruidl201ll ) to 3-D hydrodynam- 
ical simulations. Nevertheless, despite the limitations, our 
2-D models are adequate for understanding how the line 
profiles of t/a are modified by the carving of its wind by rjB- 



11 CONCLUDING REMARKS: 

IMPLICATIONS OF THE WWC CAVITY 

We presented the first 2-D radiative transfer modeling of 
the central source in Eta Carinae, taking into account the 
carving of the stellar wind of the primary star that occurs 
due to the wind of a close companion star. We investigated 
the effects of the cavity using multi-wavelength diagnostics 
around apastron of the orbital cycle. We believe that the 

2- D models presented here pave the way for future work on 

3- D, full line blanketed radiative transfer modeling of the 
Eta Car binary system, that are warranted to gain further 
insights on the effects of t)b on t/a- 

The 2-D models have several important implications for 
the understanding of the Eta Car binary system. Our mod- 
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els predict that the WWC cavity induces latitudinal and 
azimuthal dependence in the continuum and line profiles, 
which could be tested observationally with spectra reflected 
off the ejecta that surrounds the binary system. Such an in- 
vestigation might provide clues on whether or not tja is a 
rapid rotator. 

Our 2-D models unambiguously show that lines formed 
over a large volume of the wind of t/a, such as Ha, H/3, and 
Fe n lines, are affected by the cavity even during apastron. 
In addition, ionization from t)b may also affect the strength 
and shape of these lines. Therefore, analyses of long-term 
variations of Eta Car, possibly caused by changes in any 
stellar or wind parameter of one (or both) stars, should take 
into account the effects discussed in this paper. 

In addition, the presence of a WWC cavity and WWC 
zone with an orientation an d structure predicte d by hydro- 
dynamical simulations (e. g. iMadura et al.|[2012l ) causes the 
line-of-sight to the central source to have peculiar conditions. 
The combination of a half-opening angle of the WWC cavity 
~ 50° —60° and orbital inclination angle ~ 140° implies that 
the line-of-sight to the central source might cross a signifi- 
cant fraction of the post-shocked, dense regions of the WWC 
zone. As we discussed in this paper, emission and absorption 
of atomic species may occur. In addition, if dust is able to 
form along th e WWC zone, as suggested from studies of WR 
binaries (e. g. I Williams et al.ll2009l. l201ll) and Eta Car itself 
llFalceta-Goncalves et ai]|2005l ; iKashi fe SokerlboOSl ; ISmithl 
2010), emission from the central source would suffer from 
dust extinction. This scenario might explain the anoma- 
lous extinction seen tow ards the central sourc e compared to 
ejecta just ~ 0'.'3 away ( Hillier fe Alien! 1 19921 ), perhaps ex- 
plain ing the "coronagraph" effect discussed by |H06l l|Smithl 
|2010| ). Moreover, Eta Car has displayed a puzzling p hoto- 
metric brightening in the last century (e. g. lFrewH2004T ) . and 
long-term changes in the opening angle of the WWC cavity 
and/or dust formation rate along the WWC zone might be 
one of the culprits. 

Finally, we believe the 2-D models presented here pro- 
vide the foundation to understand the behavior of spectral 
lines at phases around periastron, even if 3-D models are 
ultimately needed because of the wrapping of the WWC 
cavity around r/A- Our models predict that the amount of 
P-Cygni absorption seen in H and Fe 11 lines depends on 
the amount of wind of t^aui the line-of-sight to the observer. 
Since 3-D hydrodynamical simulations show an increase in 
mate rial from r?AS wind in the line-of-sight just after perias- 
tron ( Okazaki et al 1 l2008l ; iParkin et alll2009l ; IMadura et al.l 
|2012| ). we suggest that the increase in the P-Cygni ab- 
sorption of H and F e 11 lines observed at these phases 
ilDamineh et al.lll998l ; iNielsen et ai1l2007l ; iRichardson et all 
|2010| ) is caused by the higher optical depth of these lines due 
to the increased amount of j;a's wind in the line-of-sight to 
VA- 

Below we summarize our main findings and conclusions. 

1. We made a detailed comparison of observed HST/STIS 
spectra near apastron (0 sp = 10.410 and 10.603) with 1- 
D spherical CMFGEN radiative transfer models. We find 
that a revised wind terminal velocity of «oo = 420 kms -1 
provides better agreement with the redshifted part of H 
and Fe 11 emission lines in the optical spectrum, and the 
P-Cygni absorption component of low-ionization UV reso- 



nance lines of C 11, Si 11, Fe 11, and Mg 11. We find the 
following parameters for our revised spherical CMFGEN 
model of r) A - = 420 kms" 1 , M = 8.5 x 10" 4 M©yr _1 , 
/ = 0.1, = 35200 K (at r Ross = 130), T eff = 9400 K 
(at tross = 2/3), I* = 5 x 10 6 Lq, assuming a distance of 
d — 2.3 kpc. 

2. We adapted the 2-D radiative transfer code of IBH05I to 
account for the modification of the wind of 7]a caused by the 
presence of a low-density WWC cavity and a dense WWC 
zone. We investigated in detail the effects of the WWC cav- 
ity on spectral lines, showing that two factors dominate the 
influence: the size of the line formation region compared to 
the distance of the apex of the WWC cavity, and the ori- 
entation of the cavity with respect to the observer. For our 
preferred orbital orientation with i = 138° and uj — 260°, 
the WWC cavity affects more the blueshifted and zero ve- 
locity parts of spectral line profiles. 

3. According to our 2-D radiative transfer modeling, Ha, 
H/3, and Fe 11 lines are the most affected by the WWC cavity, 
since they form over a large volume of ?;a's stellar wind. The 
excess P-Cygni absorption seen in Ha, H/3, and optical Fe 11 
in spherical models becomes much weaker or absent in 2-D 
cavity models with i — 138° and ui — 260°, in agreement 
with the HST/STIS observations around apastron. 

4. The same 2-D model simultaneously provides a superb 
fit to the observed UV spectrum of Eta Car at apastron. 
The better fits are because the UV spectrum of Eta Car is 
strongly dominated by absorption of Fe 11 lines that become 
significantly weaker when the presence of the low-density 
WWC cavity is taken into account. 

5. We obtain that orbital orientations with 110° < i < 
140° and 210° < ui < 330° are needed to simultaneously fit 
the UV and optical spectrum of Eta Car around apastron, 
assuming a in the range 50° — 70° . Larger values of a and/or 
photoionization of the wind of rjA and WWC zone by r/B 
could change the range of allowed i and uu slightly. 

6. Our models predict that the WWC cavity induces latitu- 
dinal and azimuthal dependence in the continuum and line 
profiles. In particular, the P-Cygni absorption is stronger 
in spectral lines viewed from i — 180° than at the cen- 
tral star position. Therefore, our models predict that spectra 
scattered off the Homunculus poles should have stronger P- 
Cygni absorption than those observed on the central star or 
equatorial positions. The presence of the WWC cavity thus 
mimics the effects that could, in principle, also arise in a 
single, fast-rotating t)a- 

7. Our 2-D model indicates that the apex of the WWC 
cavity does not reach optically-thick regions in the con- 
tinuum in the UV, optical, and infrared during apastron, 
and that these regions will be affected by a bore-hole ef- 
fect only at 0.95 < or b 5, 1-05. Th i s is in agreement wit h 
the findings fr o m IMadura fe Owockl <|2010h ; lMadural (|201Cj ): 
IMadura et ail l|201ll ). 

8. We find that for continuum wavelengths with A < 912 A 
or ^ ^ 80 \im, the WWC cavity penetrates the t — 2/3 
surface at apastron. For A < 912 A a pronounced bore- 
hole effect occurs, increasing by 5 to 10 orders of magnitude 
the number of photons escaping through the WWC cavity. 
However, this extra flux is easily dwarfed by that expected 
from the hot, luminous companion r/B during apastron. 

9. For 912 < A < 2000 A, even if the WWC cavity does not 
penetrate the r = 2/3 surface at apastron, a flux increase of 
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the order of 10-30% is seen in the 2-D model with a WWC 
cavity open toward the observer compared to the spherical 1- 
D model. We found that the flux increase is not only caused 
by the modification of the optical depth scale of the primary 
wind, but also by the steeper optical depth dependence of 
the source function in the UV compared to the optical/near- 
IR. 

10. We present predictions of continuum images at apas- 
tron in the UV, optical, and near-IR, that could be tested 
observationally with interferometric techniques. In addition 
to the bore-hole effect, emission and absorption from the 
dense walls of the WWC zone might significantly affect the 
continuum images and flux, depending on the orbital orien- 
tation and density of the walls. 
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